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ABSTRACT 

The cold dark matter (CDM) model of structure formation, normalized on large scales, 
leads to excessive pairwise velocity dispersions on small scales. In an attempt to circum- 
vent this problem, we study three scenarios (all with = 1) which have more large-scale 
power and less small-scale power than the CDM model: 1) an admixture of cold and hot 
dark matter; 2) cold dark matter with a non-scale-invariant, power-law primordial power 
spectrum; and 3) cold dark matter with coupling of dark matter to a long-range vector 
field. Despite the reduced small-scale power, when such models are evolved in the non- 
linear regime to large amplitude, the velocities on small scales are actually increased over 
CDM with the same value of erg. This 'flip-over', in disagreement with the expectation 
from linear perturbation theory, arises from the nonlinear coupling of the extra power on 
large scales with shorter wavelengths. However, due to the extra large-scale power, the 
recent COBE DMR results indicate smaller amplitudes for these models, ag ~ 0.5 — 0.7, 
than for CDM (for which erg ~ 1.2). Therefore, when normalized to COBE on large 
scales, such models do lead to reduced velocities on small scales and they produce fewer 
halos compared with CDM. Quantitatively it seems, however, that models that produce 
sufficiently low small-scale velocities fail to produce an adequate distribution of halos. 
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1. INTRODUCTION 



Recent observations of fluctuations in the cosmic microwave background by the 
COBE sateUite (Smoot et al. 1992) have placed the gravitational instability of dark 
matter theory of structure formation on firmer footing. At the same time, the shape of 
the linear density fluctuation spectrum is now constrained by a variety of observations 
ranging from galaxy scales (~ Mpc) up to the very large scales (~ 1000 Mpc) probed by 
COBE. The most popular theory of galaxy formation, the cold dark matter model with 
a scale-invariant spectrum of primordial fluctuations (hereafter CDM), has been studied 
by many authors using numerical A^-body experiments over a wide range of scales (e.g. 
Davis et al. 1985, hereafter DEFW; Melott 1990; Park 1990; Couchman & Carlberg 
1992; Gelb 1992, hereafter G92). Using the COBE DMR measurements to flx the density 
fluctuation amplitude on large scales, it is now necessary to re-evaluate the predictions of 
the CDM model on smaller scales. In fact, it is difficult to reconcile the COBE measure- 
ments, which indicate a rather high amplitude for CDM, with the observed low galaxy 
pairwise velocity dispersions on scales of a few Mpc. In a search for solutions to this 
problem, we consider three scenarios, each of which might be considered a variation on 
the CDM model and has less small-scale and more large-scale power than CDM. 

Given current theoretical uncertainties, the amplitudes of the initial power spec- 
tra for the models are free parameters; they are related, for example, to the self- 
coupling of the scalar field that drives infiation. It is common to normalize spec- 
tra by setting the rms density fluctuation in spheres of radii 16 Mpc to ag, assum- 
ing linear growth of modes independent of their wavenumber k = 27r/A. (We assume 
Hq = SOkms"-*^ sec~^ and 1] = 1 throughout.) Thus, for a power spectrum P{k), 



The pairwise velocity dispersion of galaxies on small scales was determined from 




{sinkR-kRcoskR) ; = 16 Mpc . (1) 







On small scales, ug is constrained by the pairwise velocity dispersion, a\\ , deflned as the 

1/2 

rms velocity along lines of separation of galaxy pairs: cr||(r) = {{v\\ — (f||))^) , where 
— {v2 — vi) ■ r/r is the parallel, peculiar velocity and f = f2 — ri- 
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redshift surveys by Davis & Peebles (1983): a\\ (r < 3 Mpc) ~ 300 ± SOkms"^ (Bean 
et al. 1983 found 250 ± 50kms~^). Researchers have compared this value with N- 
body simulations in order to constrain the CDM model amplitude a%. DEFW found 
(78 = 0.4 and G92 found ag ^ 0.5 (allowing for galaxy dispersions in simulated clusters to 
be less than virial estimates). However, the recent COBE measurements imply a larger 
amplitude for the CDM model: the observed rms fluctuation on 10°, [a-T(10°)]DMR = 
(1.085 ± 0.183) X 10-^, sets ~ 1.17 ± 0.23 [Adams et al. 1992; a nearly identical 
range for erg is obtained by fitting the COBE angular correlation function C{6)]. A larger 
amplitude (erg — l-29[li;g5], Adams et al. 1992) is also required for CDM by the large- 
scale peculiar velocities of galaxies, determined by Bertschinger et al. (1990) using the 
POTENT reconstruction algorithm. Thus, for CDM, large-scale observations suggest a 
a"8 amplitude roughly twice as large as that indicated by a"|| on small scales. 

A possible way to reconcile large scales with small scales was suggested by Couchman 
& Carlberg (1992), who found that cry for the resolved halos is a factor ~ 2 lower than 
(T|| for the mass, an effect known as 'velocity bias'. However, G92 showed that despite 
a velocity bias of = cr|| (halos)/(7|| (mass) ~ 1/2, high amplitude CDM still produces 
cr|| on scales of r > 1 Mpc in excess of the observed cr|| . Furthermore, CDM simulations 
result in an overproduction of halos, with the high-mass halos (from mergers) possibly 
representing clusters. If so, then dividing up the clusters into halos eliminates the velocity 
bias altogether (see G92). The alternative approach of using peak particles to represent 
galaxies (Bardeen et al. 1986) also leads to 6„ ^i; 1 (DEFW; Quinn, Katz, & Gelb 1992). 
We therefore claim, that ay (mass) adequately reflects (jy (halos), and consequently only 
focus on cry (mass) in the following analysis. 

2. MODELS 

We present results from 128^ particle simulations using the particle- mesh technique 
on a 256^ grid (Bertschinger & Gelb 1991) in boxes of comoving length 200 Mpc on a side. 
This size is large enough to encompass the waves that contribute significantly to and 
cry, although these simulations, with particle mass rup = 2.6 x 1O^^M0, are inadequate 
for resolving individual halos. (We use 51.2 Mpc boxes, rUp = 4.4 x 1O^M0, at the end 
of §4 to study halo formation.) 
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We explore three models with less small-scale power and more large-scale power than 
the CDM model: (1) models with an admixture of hot and cold dark matter; (2) a cold 
dark matter model with a non-scale- invariant, power-law primordial spectrum, P{k) oc 
k^, with n = 0.7; and (3) cold dark matter models in which the dark matter couples to 
a hypothetical long-range field (Frieman & Gradwohl 1991; Gradwohl & Frieman 1992). 
In all cases, we assume negligible baryon density, and, with the exception of (2), a scale- 
invariant (n = 1) primordial spectrum. In case (3), the force law between objects of mass 
m is: 

r^F/iGn?) = 1 + a(l + r/Ao) exp(-r/Ao) , (2) 

so that gravity is effectively retarded for a < in the case of a vector field. Here, a is a 
measure of the relative strength of the new force, and Aq is its range. For these 'alpha' 
models, we regenerate the optimal Green's function (Hockney & Eastwood 1982) for the 
potential corresponding to the force law of eqn. 2, 

- 1 a 

A = Xotto/a is the comoving scale at expansion factor a, and S is the Fourier transform 
of the density contrast. (We define a = 1 when erg = 1 throughout. The comoving scale 
equals the physical scale at present day expansion factor oq = erg-) Because of limited 
force resolution, the retardation effect is washed out at later expansion times — the 
principle effect enters in the initial conditions (Gradwohl & Frieman 1992) and the early 
evolution. 

Other authors have studied some of these models. Gen et al. (1992) simulated a 
erg = 0.5 n = 0.7 cold dark matter model with comparable resolution. The authors 
argued that this model has lower cr|| ~ 400 — 500 km s"-*^ than erg = 1 GDM, but that a 
velocity bias by ~ 1/1.5 is still needed to match the observed cry ~ 300 kms"-*^. However, 
we do not allow a velocity bias factor for reasons discussed earlier. Davis, Summers, 
& Schlegel (1992) performed high resolution simulations with cold and hot dark matter 
(^HDM = 0.30) in 14 Mpc boxes at as = 0.9. We demonstrate that not only is this box 
too small to adequately measure cry (as pointed out by the authors), but that large-scale 
waves (not present in a 14 Mpc box) at large values of erg can actually increase (7|| on 
small scales relative to GDM. 



(3) 
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3. LINEAR PERTURBATIONS 



We present linear realizations of the initial power spectra and linear estimates of 



are plotted in fig. la as a function of comoving wavenumber k, where Ak = (27r/L)^ 
with L = 200 Mpc. We begin all simulations at an expansion factor a = 1/50. We 
use the transfer function of DEFW for models other than the cold plus hot models, and 
the density- weighted transfer function of van Dalen & Schaefer (1991) for the cold plus 
hot models. All of the simulations use the same initial random numbers scaled to the 
appropriate power spectra. Results are shown for six simulations: CDM (standard cold 
dark matter); C+H13 and C+H30 (cold dark matter mixed with 13% and 30% hot dark 
matter); TILT7 (cold dark matter with n = 0.7); and ALPHAS and ALPHA5 (cold dark 
matter alpha models with a = —0.3, Aq = 100 kpc; and a = —0.5, Aq = 500 kpc). The 
sharp cutoff corresponds to the three-dimensional Nyquist wavenumber and the solid line 
is a direct plot of the DEFW CDM power spectrum for discrete kL /{27r). The first bin 
for the simulated models is an average over 18 waves. 

The TILT7, C+H13, and ALPHAS power spectra are similar. However, the alpha 
simulations will lag behind the others as the modified force law continues to retard the 
subsequent evolution. (We discuss the nonlinear power spectra of fig. lb in §4.) 

Although it is a poor approximation for the small scales of interest, it is nevertheless 
instructive to estimate a\\ in linear perturbation theory where the velocity and density 
are related by (Peebles 1980; Lightman & Schechter 1990) 



where 5 is the density contrast; therefore, v\d?k oc HQP{k)/k'^d^k. As a result, the linear 
estimate is 



cry for the models. The dimensionless linear power spectra P{k){Ak)^ (all with ug = 1) 



(4) 
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The factor 1 — sia{kr)/{kr) filters out the contribution of long waves to the small-scale 
pairwise velocity dispersion; for waves with A » r, the bulk flow does not contribute to 
the dispersion. This is opposite to the top hat filter Wth in eqn. 1, which filters out 
contributions from large kr. 

In fig. 2a we plot linear theory estimates of for the various spectra at ug = 1. The 
Davis & Peebles (1983) estimates from the observations are ~ 300 ± 50kms~^ for scales 
~ 1 — 3 Mpc. The models with reduced small-scale power appear to fare well at these 
scales, but, as we now show, nonlinear effects radically alter a"||. 

4. NONLINEAR CALCULATIONS 

4.1 PAIRWISE VELOCITY DISPERSIONS 

We first study nonlinear power spectra and a\\ for the models. In fig. lb we present 
nonlinear power spectra at us = 1 for several models, computed from our simulations 
in the 200 Mpc box. For low values of k, the spectra agree with their counterparts in 
fig. la, indicating a sufficiently large box size. The models with less small-scale power in 
the initial conditions continue to have less small-scale power in the nonlinear regime, yet 
the nonlinear spectral shapes differ from the linear regime and the differences among the 
various models are significantly reduced, compared to our linear estimates. 

The nonlinear a\\ for the models are shown in figs. 2b and 2c. Comparing with 
fig. 2a, it is clear that linear theory is a very poor estimator of both the amplitude and 
the general characteristics of a\\. In fig. 2b a\\ is plotted at erg =0.25, 0.5, 0.75, and 1 for 
CDM and C+H30. At low amplitude, as ^ 0.5, cry for C+H30 is lower than a\\ for CDM, 
in agreement with the expectation from linear theory. However, when the models are 
evolved further, i.e., for as ^ 0.7, a\\ 'fiips over': despite its reduced small-scale power, 
the C+H30 model yields larger cr|| than CDM. Fig. 2c shows that this fiip-over at high 
as is generic for these models; it is a reflection of the fact that their extra power on large 
scales couples significantly to small scales. This effect is manifest in the pairwise velocity 
dispersion, but not the power spectrum, because the former is more sensitive to long 
wavelengths. 

The flip-over in the initial power spectra occurs on scales exceeding ~ 75 Mpc (fig. la) 
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and is therefore missed in the 14 Mpc simulations of Davis, Summers, & Schlegel (1992); 
at as = 0.9 (the case studied there) C+H30 actually yields higher (7|| than the CDM 
model. (Originally we used 51.2 Mpc boxes and did not see this effect either.) Cen et al. 
(1992) may not have elucidated this effect since they only evolved their n = 0.7 simulation 
to as = 0.5. Our TILT7 results agree with Cen et al. (1992) at this amplitude. 

To compare these results with the observations of cr|| ~ 300 ± 50kms~-'^, we fix the 
as amplitudes of the models using the DMR data and linear perturbation theory (which 
is valid on the large scales probed by COBE). As noted in §1, for CDM, COBE yields 
(Tg ^ 1, which implies cry in excess of 1200 km s~^ over scales of a few Mpc. For the other 
models, COBE's rms ffuctuations on 10°, including the DMR errors and cosmic variance 
for the models, imply as ^ 0.78 ±0.16 (for C+H13), 0.69 ±0.14 (for C±H30), 0.53 ±0.11 
(for TILT7), 0.95 ±0.19 (for ALPHA3), and 0.51 ±0.10 (for ALPHA5). When normalized 
to COBE, it is possible to find models with reduced cry. One should, however, point out 
that models with cr|| ~ 400 — 550kms~^, although favorable over erg = 1 CDM, may still 
be inadequate. This is due to the fact that simulation-to-simulation variations in ay are 
typically ^ lOOkms"^ (G92), the observed errors in (jy are ~ 50kms~^, and Bean et al. 
(1992) found lower cry of order 250 ± 50kms~^ than Davis and Peebles (1993). 

In fig. 3a we plot cry for 'favorable' models, subject to the above COBE normaliza- 
tions. They are: erg = 0.6 ALPHA5 and ag = 0.5 TILT7, both consistent with COBE; 
as — 0.5 C+H30, which is somewhat beyond the la level; and as — 0.5 CDM, which is 
inconsistent with COBE and only shown here for comparison. We exclude C+H13 and 
ALPHA3 from our list, as they produce excessive (7||. TILT7 and C±H30 both produce 
small-scale cry in the 400 — 550kms~^ range, a factor of ~ 2 smaller than erg = 1 CDM, 
but still at least 30% above the observed uy. The only model in fig. 3a that matches 
COBE and small-scale cry is ALPHA5. It nicely reproduces the observed cry, but as we 
now demonstrate, suffers from inadequate halo formation. 

4.2 HALO FORMATION 

In addition to the affect on cry, reducing small-scale power can also help alleviate 
some problems associated with excessive low-mass halo formation in the field and an 
excessive number density of high-mass objects that plague the CDM model (White et 
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al. 1987; G92). In order to study halo formation, we simulate models in 51.2 Mpc 
boxes, again using 128^ particles and a 256^ particle-mesh grid. (We emphasize again, 
that although adequate for analyzing the halo distribution, the 51.2 Mpc box simulations 
cannot be used to compute cr||.) 

The details of galaxy formation require higher resolution simulations with separate 
hot and cold dark matter particles (e.g. Davis, Summers, & Schlegel 1992) and separate 
dark matter and baryonic matter particles in alpha models (the coupling to the vector 
field only occurs for the dark matter and can lead to a natural bias between dark and 
baryonic matter, Gradwohl & Frieman 1992). These considerations, however, are not 
likely to significantly affect mass pairwise velocity dispersions. 

Halos in our simulations are identified, by use of the DENMAX algorithm 
(Bertschinger & Gelb 1991), as local density maxima in the evolved, nonlinear density 
field. The distribution of halos, characterized by their circular velocities (computed from 
the enclosed mass within R =300 kpc of the DENMAX center, i.e. ^/GM{<K)JR) are 
then counted in 25kms"^ bins. The results are shown in fig. 3b for the same scenarios 
as in fig. 3a. The solid line is the observed estimate, using a Schechter (1976) luminosity 
function coupled with TuUy-Fisher (1977) and Faber-Jackson (1976) relationships (as- 
suming that 70% of the halos are spirals and 30% are ellipticals) to relate luminosity to 
mass. (See G92 and White et al. 1987 for details. The observed estimates as shown tend 
to overestimate the highest-mass halos compared with complete elliptical surveys, and 
the overall normalization has an error ~ 30%.) 

The cases erg = 0.5 CDM and = 0.5 TILT7 appear to match the observed distri- 
bution fairly well. G92 demonstrated, however, that high-mass halos should be divided 
into clusters of halos so that 1) the simulations contain clusters and 2) extra weight is 
given to dense systems thereby enhancing the two-point correlation function in biased 
(cTg < 1) models (White et al. 1987). Therefore, TILT7 actually does much better than 
CDM — it produces less mid-mass halos and less high-mass halos which, when divided 
into mid-mass cluster members, make up the mid-mass deficit. C-I-H30 produces too few 
halos, which can be remedied by evolving the simulation further at the expense of raising 
cry. ALPHA5 drastically fails to match the observed distribution for erg = 0.6. (ALPHA3 
does better, but it produces excessive cr||.) 
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5. CONCLUSIONS 

For as ^ 0.5 (the precise value depends on the model), the nonlinear coupling of 
waves in models with more large-scale power and less small-scale power than CDM ac- 
tually increases (7|| on small scales, in complete disagreement with linear perturbation 
theory. At ag = 0.5, CDM, C+H30, and TILT7 yield a\\ ~ 400 - 550 kms-^ on smaU 
scales, but CDM has inadequate large-scale amplitude. C-I-H30 and TILT7 also gener- 
ate, consistent with observations, less halos than CDM. The only model which produces 
cr|| ~ 300kms~^, and at the same time matches COBE, is ALPHA5 at ag — 0.6, but it 
fails to produce a sufficient number density of halos. 

The fact that it seems difficult (if not impossible?) to accommodate a low , and 
still have enough small-scale power for adequate halo formation, may hint to a basic 
problem of Q = 1 cosmogonies. One way out of this apparent impasse is, of course, to 
lower the matter density of the universe (recall that v oc O"^/^), and thereby maintain a 
low cr|| with increased small scale power. It is clearly too early to view this problem as a 
death stroke to = 1 scenarios, more work still needs to be done involving the details of 
halo formation and biasing. 

We summarize our results in Table 1. There are no obvious 'winners'. The balance 
of scores can be shifted by varying as', e.g., C+H30 can be shifted to +, — ,+ by choosing a 
higher amplitude. Recent COBE measurements on large scales require the investigation 
of a myriad of models with more free parameters than CDM. On small scales, in = 1 
scenarios, it seems difficult to reconcile low a\\ with sufficient halo formation, and one 
may be forced to consider scenarios with f2 < 1 or a cosmological constant. 
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Corporate Research Institute. This work was supported in part by the DOE and NASA 
grant NAGW-2381 at Fermilab. B.C. was also supported by DOE #DE-FG03-91ER 
(40662 Task C) at UCLA. 
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TABLE 1 
Scorecard 



Model 


(78 


COBE 


C7||(r < 3 Mpc) 


Halos 


CDM 


0.5 




o 




C+H30 


0.5 


+ 


o 


o 


ALPHAS 


0.6 


+ 


+ 




TILT7 


0.5 


+ 


o 


+ 



Table 1 Caption: Scorecard for the various models with + indicating a favorable score, 
o indicating a marginal score, and — indicating a disfavorable score. 
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FIGURE CAPTIONS: 



Figure 1: a) Linear realizations (all use equivalent random numbers) of various power 
spectra normalized to cts = 1. (C+H30 has less small-scale power than C-I-H13 and 
ALPHA5 has less small-scale power than ALPHAS.) b) Nonlinear power spectra at 
(Tg = 1 for some of the simulations shown in fig. la. In both a) and b), the solid curve is 
the analytic form of the DEFW linear CDM power spectrum. 

Figure 2: a) Linear theory calculations of (T|| versus comoving separation (eqn. 5) for 
as = 1. (C+H30 is lower than C+H13 and ALPHAS is lower than ALPHA3.) b) 
Nonlinear calculations of cr|| for CDM (short dashed curves) and for C-I-H30 (long dashed 
curves) at four values of cxg. c) Nonlinear calculations of a\\ for CDM (short dashed 
curves), C+H13 (long dashed curves), ALPHA3 (dot-short dashed curves), and TILT7 
(dot- long dashed curves) at (T|| = 0.5 and 1. The observed estimates from galaxies for 
r ~ 1 -3 Mpc are SOOiSOkms"^ (Davis & Peebles 1983) and 250±50kms-^ (Bean et 
al. 1983). 

Figure 3: a) Nonlinear calculations (200 Mpc box) of cry. b) Distribution of halos versus 
circular velocity (in AV^irc = 25kms~^ bins). The solid curve is the observed estimate. 
In both a) and b) the cases are: CDM at ag = 0.5 (short dashed curves), C-I-H30 at 
cTg = 0.5 (long dashed curves), ALPHA5 at ag = 0.6 (dot-short dashed curves), and 
TILT7 at as = 0.5 (dot- long dashed curves). 
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